Extended Hubbard model of superconductivity driven by charge fluctuations in 

iron-pnictides 
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We present a scenario for iron-pnictide superconductivity mediated by charge fluctuations tfiat 
are strongly enhanced by Fe-As intersite electronic interactions. Deriving an eight-band extended 
Hubbard model including Fe 3d and As 4p orbitals for the LaOFeAs family, we show that charge fluc- 
tuations induced by p-d charge transfer and As orbital polarization interactions in the Fe-pnictogen 
structure peak at wavevectors (0,0), and (tt, 0) and (tt, tt) respectively. Intraorbital spin-singlet 
pairing attraction develops at these wavevectors and the solution of the linearized gap equation 
shows robust s-wave superconductivity with both s± and s++ gap functions. 

PACS numbers: 74.70.Xa, 74.20.Mn, 74.20.Rp, 74.20.-z 



The mechanism of high-Tc superconductivity in the 
Fe-pnictides has attracted enormous attention since its 
original discovery in F-doped LaFeAsO (1111) (H- The 
majority of the theoretical efforts has focused on the 
proximity of the superconducting (SC) phase to the spin 
density wave (SDW) state and the multiple Fermi sur- 
faces (FS) associated with the Fe 3d and As Ap orbitals 
0-01 ■ emerging picture is that spin fluctuations and 
FS scattering favor spin-singlet s±-wave pairing where 
the gap function changes sign between hole and elec- 
tron FS due to the intraorbital repulsion in the particle- 
particle channel. For the prototypical 1111 series, where 
Tc reaches the record high of 55K when La is replaced by 
other rare earths [1] , NMR Knight shift measurements in- 
deed find spin-singlet pairing but it remains unclear 
whether spin fluctuations are the driving force for super- 
conductivity. Upon electron-doping, spin fiuctuations in 
the normal state are dramatically suppressed; the SDW 
phase terminates abruptly and is separated from the SC 
state by a first order-like transition l3-l3- The corre- 
lation between Tc and the low energy spin-fluctuations 
measured by the spin-lattice relaxation rate has been 
found to be rather weak l3"12|. Moreover, applying 
pressure near the optimal doping level increases from 
23K to 43K while the strength of spin-fluctuations re- 
mains unchanged fl^ . This is further supported by re- 
cent muon spin rotation (^SR) and magnetization exper- 
iments in the underdoped regime that show hydrostatic 
pressure suppresses magnetic interactions but strongly 
enhances Tc [l5| . Thus, spin- fluctuations alone cannot 
fully account for the pairing mechanism of iron-pnictide 
superconductors . 

In this paper, we explore a different scenario where 
the superconductivity is driven by charge fluctuations. 
There are indeed emerging experimental evidence that 
the pnictides are close to the charge ordering instabil- 
ity. In the 1111 series, two distinct charge environments 
are detected by As NQR measurements in the under- 
doped regime, indicative of local electronic charge order 



3 12 1- In contrast to the cuprates, the Fe-pnictides are 
p-d charge transfer metals with low energy charge fluc- 
tuations. It is thus important to go beyond the local 
Hubbard interactions and consider the interatomic inter- 
actions. Furthermore, due to the large spatial extent of 
the As 4p orbital, the interactions between the Fe 3d and 
As Ap electrons are important both in the charge trans- 
fer channel and in the As orbital polarization channel 
when charges fluctuate at the Fe site. We found that 
it is a generic feature of the Fe-pnictogen structure that 
these interactions produce enhanced charge fluctuations 
at (0,0), (7r,7r), and (tt, 0) respectively, and mediate at- 
tractions for intraorbital pairing at these wavevectors. 

We focus on the electron-doped 1111 series that shares 
a single Fe-pnictogen layer per unit cell and is the most 
quasi-two-dimensional Fe-pnictides. We derive an ex- 
tended Hubbard model as the low energy effective Hamil- 
tonian for the FeAs layer: H = Hq + Udd + Vpd, where 
Hq is a tight-binding model for the band structure in- 
cluding both the Fe 3d and As ip orbitals; Udd describes 
the local interactions, intraorbital Hubbard repulsion U 
and Hund's rule coupling J, at the Fe sites; and Vpd con- 
tains the nearest neighbor (NN) charge transfer interac- 
tion V and As orbital polarization interactions AVi for 
Px-Py and AV2 for Pz-Px,y Treating U and J as effective 
interaction parameters, a random phase approximation 
(RPA) study of the charge and spin fluctuations is car- 
ried out as a function of the Fe-As inter-site interactions. 
We find that the enhanced charge fluctuations lead to 
robust s-wave superconductivity with both s± and 5-1-4. 
gap symmetry as summarized in Table I for a wide range 
of doping levels. 

The low energy part of the Lallll band dispersions 
shown in Fig. la can be described by a tight-binding 
model Hq for the Fe 3c? and As Ap complex [7|. For the 
single-layered 1111, it is possible to unfold the reduced 
zone to the original one corresponding to one FeAs per 
unit cell and work with 8 bands specified by an orbital 
index a = l{dxy), '^{dyz), 3(4^;), 4{d^2_y2), 5(d3^2_r2), 
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TABLE I: Symmetry of the leading pairing instability driven 
by p-d interactions for different on-site U and Hund's rule 
coupling J. All cases listed are nodeless. 
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FIG. 1: (Color online) Eight-band p-d model (a) Comparison 
of the band dispersions to LDA band structure in the reduced 
zone. Line thickness and symbol size denote Fe 3d content. Fe 
3d (b) and As 4p (c) contributions to the FS in the unfolded 
zone at 10% electron doping. Symbol sizes denote the orbital 
content with those of the As 4p enhanced by a factor of 4. 

6(px), 7(py), 8{pz)- Fig. la shows that the p-d model 
Hq describes well both the LDA band dispersion and the 
orbital character for the undoped case with 12 electrons 
per unit cell. At 10% electron doping, the FS contain 
two hole pockets (labeled by a and /3) centered around 
r and two electron pockets around X (labeled by 7) and 
Y pionts. Fig. lb and Ic display the dominant Fe 3d and 
As 4p orbital characters on the FS respectively. 
The electronic interactions have the general form 

ij,aa' ab^a'b' 

where cl^^ creates a spin-cr electron on site i in orbital a. 
The Coulomb integral is given by 

Wab,a'b'{r^J)^ I dhd^'v' cj):{r)<j>l, {r')V {R,j)^a' {r')Mr) , 



(2) 

where Rij = jry + r' — r| and (jia is the wavefunction 
of orbital a. Retaining the dominant on-site interactions 



for the Fe atoms (those of the As are much weaker) and 
the NN p-d interactions, we write Hj = Udd + Vpd- Udd 
attains the usual multi-orbital Hubbard model 

Udd^U'^nia'fnial+ (U' -^Jj fltafliP (3) 

i.a ^ ^ i,a<(i 

with intra and inter orbital on-site Coulomb repulsions 
U — Waa,aa{0), U' — Waa,iip{^) and the Huud's rule 
coupling J = J' = Wap.api^)- Orbital rotation symme- 
try requires U ~ U' -\- 2 J. Here and henceforth, we use 
a, (3 = 1, 2, • • • ,5 to distinguish Fe 3d orbitals from As 
4p orbitals denoted hy — 6, 7, 8. 

The Coulomb integral Wap,^v describes a rich variety 
of Fe-As interatomic interactions. The p-d charge trans- 



fer interaction Va 



Wo. 



(t*) where r* is the vector 



connecting the NN Fe and As. The importance of V was 
emphasized in the context of the cuprate superconduc- 
tivity [ll]. Furthermore, l^Va.^u — W^aa,/ii/(r*) describes 
the As 4p orbital polarization induced by the Fe elec- 
tric field associated with the charge fluctuations in the 
a-orbital. This is different from the hi gher energy As 
4p-5s polarizations discussed in Ref. fiol |2C| . The large 
spatial extent of the As 4p orbital has important conse- 
quences: (i) The bare interaction AVi^2 estimated using 
the hydrogen-like atomic wavefunctions in Eq. ([2]) is re- 
markably large and about 10-20% of the p-d charge trans- 
fer V . Since V is subject to charge screening whereas 
AVi,2 is not, the effective interaction strengths can be 
comparable, (ii) The interaction involving the polariza- 
tion of the smaller Fe orbitals Wa^^^^ and the interaction 
between the Fe and As polarization clouds WajS^^v are at 
least one or two orders of magnitude smaller and can 
thus be neglected, (iii) Since the 3d orbitals are much 
smaller, their dependence in V and ISV can be ignored. 
We thus arrive at the following Hamiltonian for the p-d 
interactions. 



A1/2 



T- ■ n 



h.c. 



(4) 



where nf and are the total density operators of the d 
and p electrons respectively. Since the FeAs block devi- 
ates from the ideal tetrahedron structure, two interaction 
parameters, AVi and AV2, are introduced to distinguish 
between As Px-Py and Pz-Px,y orbital polarizations. Note 
that the polarization (quadrupole) term is orientation- 
dependent and T^j" accounts for the sign of the wave- 
function overlap, 
tion reads 



In momentum space, the p-d interac- 



(5) 
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2%' 

i4AV2 



where the form factors F^^(q) = cos^qx cos ^ 
■F67(q) = -iAVi sinigj; sin^qy, Fes{q) = 
sinig^, cosily, and F7g(q) = ~iAAV2 cos iq^: sinig^. 

We next present a complete RPA treatment of the in- 
teractions in Eqs.® and (jH). The charge and spin sus- 
ceptibihties can be written as 34 x 34 matrices 

r(q,^i) = X°(q,^0/[1 - U'f{q,iOi)], (6) 



where the bare susceptibihties Xab a'b'i^^^i) ~ 
-(r/iV)Ek,™G^(k + q,e™ + UjOGi.iKe^) with 
the noninteracting Green's function G'^(k, e^) — 
[iem — i?o(k)]~^. In Eq. ([5]), the nonzero elements of the 
interaction matrices t/*, U'^, and V'^ are: 



= u', u, 



C^^a,/3^ = 2C/'- J, C/^^^^„ = J', and 
The on-site interaction enhances 



(q) = F,M) 

(reduces) the spin (charge) susceptibility. The inter-site 
p-d interaction V"^, on the other hand, affects only 
the charge sector, entering in the block-off-diagonal 
elements in the denominator. They lead to enhanced 
charge fluctuations at wavevectors where the interactions 
Ffii^il) in Eq- ® a-re maximum in momentum space, 
i.e., at Q = (0,0) for V; (7r,7r) for AVi, (tt, 0) and (0,7r) 
for AV2. 

We shall describe our results for 10% electron doping 
with a moderate effective U — leV, but a reasonably 
large ratio J/U = 0.3 in accord with the large Hund's 
rule coupling in the pnictides (last row in Table I). Sev- 
eral prominent intraorbital static charge susceptibilities 
Xaa,aa(q) are shown for ^=0.26 eV (Fig. 2a), Ayi=0.3 
eV (Fig. 3a), and Ay2=0.28 eV (Fig. 4a), independently. 
Clearly, the inter-site interactions enhance the intraor- 
bital charge fluctuations by introducing peaks at the cor- 
responding Q that grow with increasing V and AVi^2- 
We verified that their emergence is tied to the soften- 
ing of the collective modes in the imaginary part of the 
dynamical charge and charge transfer susceptibility 21 1 . 
Note that the p-d interactions in Eq. ^ leave the Fe 3d 
interorbital susceptibility pa unchanged. 

To study superconductivity, we evaluate the pairing 



vertex dressed by the spin and charge fluctuations [22, 
[23I I . The effective spin-singlet pairing interaction is given 
by 

nq) =\u' + \u'r{ci)u' + ^[^' + 

Ape + 21/^(q)]r(q)[?7^ + 2y^(q)], (7) 

where x*''^(q) = x^''^{(\,uji = 0) are the static spin and 
charge susceptibilities. The spin-triplet pairing turns out 
to be sub-leading. The calculated P(q) are shown in 
Figs. 2b, 3b, and 4b for interactions V, AVi, and AV2 
respectively. Remarkably, with the enhancement of the 
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FIG. 2; (Color online). Effects of p-d charge transfer V at 
{U, J) = {1, 0.3) eV. (a) Intraorbital RPA charge susceptibility 
and (b) Singlet intraorbital pairing interaction at 1^=0.26 eV. 
(c) s- and d-wave eigenvalues A as a function of V . (d) s-wave 
gap symmetry function along three FS sheets at V^=0.264 eV 
where As = l. Angles are measured from x-axis. 

charge fluctuations near Q (peaks), the repulsion is weak- 
ened (dips) in the intraorbital pairing potential Paa,aa 
and turns into attraction for intraorbital pairing near 
Q when the corresponding p-d interaction is sufficiently 
strong. This is in contrast to the pairing interactions 
mediated by spin fluctuations that are repulsive at all 
q. The SC instability can be obtained by solving the 
linearized gap equation. 



AA,fc(k)=--^ ^ Paa.,b.fc(k-k') (8) 

k',n a'h' .a"h" 
X G°„,, (k', UJ.n)Gl„,, (-k', -LO.n)Aa.b- (k') 



in the orbital basis, where Ao6(k) is an 8 x 8 normalized 
gap symmetry function. The pairing instability sets in 
when the largest eigenvalue A reaches unity at T = Tc. 
To overcome the finite-size effects, we solved Eq. self- 
consistently at T= 20 meV on an 80x 80 momentum mesh 
to obtain A and A(j{,(k) as a function of V and AT4,2- The 
gap symmetry function can be easily transformed into the 
band basis by a unitary rotation and plotted along the 
FS. 

Superconductivity driven by inter-site interaction V 
is summarized in Fig. 2. The eigenvalues A plotted 
as a function of V in Fig. 2c show that s-wave pair- 
ing is more favorable than pairing with d-wave symme- 
tries and superconductivity sets in at a reasonably small 
Vc = 0.264eV. The normalized gap symmetry function in 
Fig. 2d shows that the pairing symmetry is the nodeless 
s-l--wave; with opposite signs for the pairing gaps on the 
electron (7) and the hole [a and /?) pockets. The ob- 
tained Aa6(k) in the orbital basis shows that all orbitals, 
including those of the As 4p, contribute in a complicated 
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FIG. 3: (Color online) Effects ofpx-py orbital fluctuation AVi 
at (U, J) — {1, 0.3) eV. (a) Intraorbital RPA charge susceptibil- 
ity and (b) Singlet intraorbital pairing interaction at AVi=0.3 
eV. (c) s- and d-wave eigenvalues A as a function of AVi. 
(d) s-wave gap symmetry function along three FS sheets at 
1^=0.309 eV where As = l. Angles are measured from a:-axis. 

manner to the behavior of the gap function on the FS. 
Nevertheless, the pairing symmetry can be quahtatively 
understood from the dominant intraorbital pairing inter- 
actions shown in Fig. 2b. While the increasing attrac- 
tion peaked around (0, 0) provides the main pairing force 
through forward scattering in contrast to spin fluctua- 
tion mediated pairing, the scattering by the repulsion 
near (tt, 0) and (0, tt) favors a sign change between the 
electron and the hole pockets in a similar manner as in 
the spin fluctuation scenario 0, [EH- Furthermore, the 
repulsion near (tt, tt) causes a degree of frustration for 
the s±-pairing, leading to the large asymmetry of the gap 
function and large variations on the electron FS. Remark- 
ably, keeping the same ratio J/U — 0.3, but reducing the 
Hubbard C/ by a factor of two, we find that the pairing 
symmetry changes to the s+^-wave due to the reduction 
in the repulsion at finite momenta associated with spin- 
fluctuations. The change from s± pairing at large U to 
s_|__)_ pairing at small U is also true for a smaller ratio of 
J/U = 0.1 and may be generic of the SC phase driven by 
the p-d charge transfer interaction V (Table I). 

Superconductivity driven by inter-site interaction AVi 
is summarized in Fig. 3. The largest eigenvalues of the 
gap equation plotted in Fig. 3c show that s-wave pairing 
dominates over d-wave symmetries and the SC phase sets 
in at AVi,c — 0.309eV. The gap symmetry function over 
the FS shown in Fig. 3d reveals that the pairing symme- 
try is the sign-changing s±-wave. Remarkably, the gap 
over the electron pocket oscillates moderately around a 
value that is close in magnitude to that on the inner hole 
pocket, but larger than that on the outer hole pocket, 
in excellent agreement with the gap ratios observed by 
ARPES in optimally doped Kj:Bai_^Fe2As2 HI- More- 



FIG. 4: (Color online). Effects of p^-px.y orbital polarization 
AV2 at {U, J) = (l, 0.3) eV. (a) Intraorbital RPA charge sus- 
ceptibility and (b) Singlet intraorbital pairing interaction at 
AV2=0.28 eV. (c) s- and d-wave eigenvalues A as a function 
of AV2- (d) s-wave gap symmetry function along three FS 
sheets at AV2=0.288 eV where As = l. Angles are measured 
from X-axis. 

over, we find that the nodeless s± pairing symmetry is 
a robust feature of the superconductivity driven by Fe 
charge fluctuations coupled to As Px-Py orbital polariza- 
tion for different values of U and J/U as shown in Table 
I. This remarkable feature is a result of the pairing inter- 
action shown in Fig. 3b. The repulsion at (tt, tt) has been 
turned into the growing attraction by AVi that provides 
the main pairing force through (tt, 7r)-scattering, leaving 
the repulsion at (tt, 0) and (0, tt) unfrustrated that locks 
the opposite sign of the gap functions on the electron and 
hole pockets. 

Superconductivity driven by inter-site interaction AV2 
is summarized in Fig. 4. It is clear from Fig. 4c that 
the leading SC instability remains in the s-wave channel 
and sets in at AV2.C = 0.288eV. The pairing interaction 
in Fig. 4b shows that AV2 has turned the repulsion at 
(tt, 0) and (0, tt) due to primarily spin-fluctuations into 
the growing attraction which serves as the dominate pair- 
ing force in this case. As a result, the s± symmetry be- 
comes unfavorable. Indeed, the gap symmetry function 
shown in Fig. 4d reveals an anisotropic s-|_4--wave with 
significant variations on the electron pocket. We find 
that the s++-wave pairing is a robust feature of the su- 
perconductivity driven by AV2 for different values of U 
and J/U, as shown in Table I. 

In summary, we proposed that the iron-pnictides su- 
perconductivity can be driven by charge fluctuations. 
The inter-site interactions in the Fe-pnictogen structure 
are found to produce strong charge fluctuations that 
mediate attractions in the spin-singlet pairing potential 
around wavevectors (0,0), (7r,7r), and (tt, 0). For elec- 
tron doped LaFeAsO, moderate Fe-As intersite inter- 
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action strengths can induce superconductivity with ro- 
bust s-wave symmetry; both sign-changing s± and sign- 
preserving S++ gap functions are possible. We suspect 
that electron-phonon couphng [25| may play a role in such 
a pairing mechanism, particularly because these wavevec- 
tors are the same as the possible lattice instability vec- 
tors. It is also tempting to speculate that the 1x2 and 
\/2 X a/2 structures observed by STM in (Ba,Sr)Fe2As2 
[26j are related to the strong As orbital fluctuations in 
the bulk pinned by the surface potential. The strong 
charge fluctuations can be pinned by impurities and de- 
fects in the bulk of the sample, leading to local charge 
order and/or orbital polarization that should be observ- 
able to local probes such as NMR and ^SR and serve as a 
test of the present theory through their correlations with 
the SC transition temperature. 
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